Suppression pattern of neutral pions at high transverse momentum in Au+Au 
collisions at v /s/v r /v=200 GeV and constraints on medium transport coefficients 
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For Au+Au collisions at 200 GeV we measure neutral pion production with good statistics for 
transverse momentum, px, up to 20GeV/c. A fivefold suppression is found, which is essentially 
constant for 5<pt<20 GeV/c. Experimental uncertainties are small enough to constrain any model- 
dependent parameterization for the transport coefficient of the medium, e.g. (q) in the parton 
quenching model. The spectral shape is similar for all collision classes, and the suppression does 
not saturate in Au+Au collisions; instead, it increases proportional to the number of participating 
nucleons, as TVpart 2 '' 3 . 

PACS numbers: 25.75.Dw 



Large transverse momentum (pr) hadrons originate 
primarily from the fragmentation of hard scattered 
quarks or gluons. In high energy p+p collisions this is 
well described in the framework of perturbative QCD [lj . 
In ultra-relativistic heavy ion collisions such hard scatter- 
ings occur in the early phase of the reaction, and the tran- 
siting partons serve as probes of the strongly interacting 
medium produced in the collisions. Lattice QCD pre- 
dicts a phase transition to a plasma of deconfined quarks 
and gluons, which induces gluon radiation from the scat- 
tered parton and depletes hadron production at high pt 
("jet quenching") [2, Q- The measurements in Au+Au 
collisions at RHIC showed suppressed hadron yields in 
central collisions 0| as predicted 0, Q , and motivated an 
advanced theoretical study of radiative energy loss using 
different approaches 0]. 

All the energy loss models must incorporate space- 
time evolution of the medium, as it is not static, as 
well as the initial distribution of the partons through- 
out the medium. Models generally also include an in- 
put parameter for the medium density and/or the cou- 
pling. Different assumptions in the various models lead 
to similar descriptions of the 7r° suppression with differ- 
ent model-dependent parameters [3, |9( . For instance, the 
Parton Quenching model (PQM) is a Monte Carlo using 
the quenching weights from BDMPS [5] that combines 
the coupling strength with the color-charge density to 
create a single transport coefficient, often referred to as 
(q) 0, El, which gives the average squared transverse 
momentum transferred from the medium to the parton 
per mean free path. 

Measurement of identified particles up to the highest 
possible pt, establishing the magnitude, pt and central- 
ity dependence of the suppression pattern, is crucial to 
constrain the theoretical models and separate contribu- 
tions of initial and final state effects from the energy loss 
mechanism. Collision centrality is related to the average 
pathlcngth of the parton in the medium. The suppres- 



sion of 7T° puts important constraints on calculations of 
the energy loss, as neutral pions can be identified up to 
very high pt- Whereas it has also been shown that the 
di- hadron suppression at high may be somewhat more 
sensitive than single hadron suppression to the medium 
opacity such improvement is contingent upon the 
theoretical and experimental, statistical and systematic 
uncertainties. 

This Letter reports on the measurement of neu- 
tral pions up to px=20GeV/c in Au+Au collisions at 
^/snn— 200 GeV at RHIC, using the high statistics col- 
lected in Run-4. Based upon the data we extract the 
(q) parameter of the PQM model for the most central 
collisions. 

The analysis used 1.03 x I0 9 minimum bias events 
taken by the PHENIX experiment [13]. Collision cen- 
trality was determined from the correlation between the 
number of charged particles detected in the Beam-Beam 
Counters (BBC, 3.0<|?7|<3.9) and the energy measured 
in the Zero Degree Calorimeters (ZDC). A Glauber model 
Monte Carlo with a simulation of the BBC and ZDC 
responses was used to estimate the associated average 
number of participating nucleons ((iVpart)) and binary 
nucleon-nucleon collisions ((JV co ii)) for each centrality 
bin Q. 

Neutral pions were measured in the ir° — + 77 decay 
channel with the photons reconstructed in the Electro- 
magnetic Calorimeter (EMCal) located in the two central 
arms of PHENIX (\r)\ < 0.35). The EMCal consists 
of two subsystems: six sectors of lead-scintillator sand- 
wich calorimeter (PbSc) and two sectors of lead-glass 
Cerenkov calorimeter (PbGl) at the radial distance of 
about 5 m. The fine segmentation of the EMCal (Sip x Srj 
~ 0.01 x 0.01 for PbSc and - 0.008 x 0.008 for PbGl) en- 
sures that the two photons from a 7r° — > 77 decay are well 
resolved up to p T » 12 (PbSc) and 16 (PbGl) GeV/c. 
Data from the two subsystems were analyzed separately 
and the fully corrected results were combined. 
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FIG. 1: Left: tt° invariant yields for different centralities (PbSc and PbGl combined). Right: consistency between the results 
obtained separately from PbSc and PbGl 



Details of the analysis including extraction of the raw 
7T° yield, correction for acceptance, detector response 
(energy resolution, dead areas), reconstruction efficiency 
(particle identification cuts) have been described else- 
where [ItJ ■ In this analysis the higher range re- 
quired correction for losses in the observed (raw) 7r°s due 
to "cluster merging" . 

With increasing 7T° momentum, the minimum open- 
ing angle of the two decay photons decreases, and even- 
tually they will be reconstructed as a single cluster. 
Such "merging" reaches 50% of the total raw yield at 
p T =14 GeV/c in the PbSc and at p T =18 GeV/c in the 
PbGl due to their different granularity and Moliradius. 
Merged clusters were rejected by various shower profile 
cuts, and the loss was determined by simulated single 7r°s 
embedded into real events and analyzed with the same 
cuts. The loss increases slowly with centrality. The sys- 
tematic uncertainties were estimated by comparing ir° 
yields in the PbSc extracted in different windows of asym- 
metry \E yi — E 12 \/(E 11 + E 12 ) and also by comparing 
yields in the PbSc and PbGl. 

We considered two sources of 7r°s not originating from 
the collision vertex: those produced in nuclear interac- 
tions of hadrons with detector material (instrumental 
background) and feed-down products from weak decay 
of higher mass hadrons (physics background). Based 
upon simulations both the instrumental background and 
feed-down background were found to be negligible (<1% 
above pT>2.0GeV/c) except for the contribution from 



TABLE I: Summary of the systematic uncertainties on the n° 
yield extracted independently with the PbSc (PbGl) electro- 
magnetic calorimeters. All but one, off- vertex 7r°, are uncor- 
rected between PbSc and PbGl, and centrality dependence is 
negligible. The last row is the total systematic uncertainty on 
the combined spectra. Detailed description of how the errors 
are correlated as a function of pt can be found in [171 ]. 



p T (GeV/c) 


2 




6 


10 




16 


uncertainty source 




PbSc (PbGl) 






yield extraction (%) 


3.0 (4.1) 


3.0 


(4.1) 


3.0 


(4.1) 


3.0 


(4.1) 


PID efficiency (%) 


3.5 (3.9) 


3.5 


(3.5) 


3.5 


(3.7) 


3.5 


(3.9) 


Energy scale (%) 


6.7 (9.0) 


8.0 


(9.2) 


8.0 


(8.2) 


8.0 


(12.3) 


Acceptance (%) 


1.5 (4.1) 


1.5 


(4.1) 


1.5 


(4.1) 


1.5 


(4.1) 


7T° merging (%) 


"(-) 




(") 


4.4 (-) 


28 


(4.8) 


Conversion (%) 


3.0 (2.5) 


3.0 


(2.5) 


3.0 


(2-5) 


3.0 


(2.5) 


off-vertex tt° (%) 


1.5 (1.5) 


1.5 


(1.5) 


1.5 


(1.5) 


1.5 


(1.5) 


Total (%) 


8.7 (12) 


9.8 


(11) 


11 


(11) 


30 


(15) 


PbSc and PbGl 
combined: Total (%) 


7.0 




7.5 




7.6 




14 



Kg decay («3% of ir° yield for pt>1 GeV/c), which has 
been subtracted from the data. Finally the yields were 
corrected to the center of the px bins using the local 
slope. 

The main sources of systematic uncertainties are yield 
extraction, efficiency corrections, and energy scale, none 
of which exhibit a significant centrality dependence. The 
PbSc and PbGl detectors have quite different systematics 
with all but one of them (off- vertex 7r°) uncorrelated. 
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FIG. 2: Nuclear modification factor (Raa) for 7r°s. Error 
bars are statistical and p-r-uncorrelated errors, boxes around 
the points indicate px-correlated errors. Single box around 
Raa = 1 on the left is the error due to JV co n, whereas the single 
box on the right is the overall normalization error of the p+p 
reference spectrum. 



Therefore, when combining their results, the total error is 
reduced in the weighted average of the two independent 
measurements. The final systematic uncertainties (one 
standard deviation) on the spectra are shown in Table [U 



: PHENIX 7t° (Au+Au^5%-eentral) 
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FIG. 3: Left: tt° Raa for the most central (0-5%) Au+Au 
collisions and PQM model calculations for different values of 
(q). Red curve: best fit. Right: x 2 (ei,,e c ,p) distribution for a 
wide range of values of (q) . 

The left panel of Figure [1] shows the ir° invariant yield 
spectra for all centralities as well as minimum bias, com- 
bined from the independent PbSc and PbGl measure- 
ments. In the overlap region the results are consistent 
with those published earlier [l6| while the errors are re- 
duced by a factor of 2 to 2.5. The right panel shows 
the ratios of PbSc and PbGl spectra to the combined 
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FIG. 4: Integrated nuclear modification factor (Raa) for ty° 
as a function of collision centrality expressed in terms of N palt . 
The error bars/bands are the same as in Fig. The last two 
points correspond to overlapping centrality bins, 0-10% and 
0-5%. The dashed lines show the fit to a function. See text. 



one at three centralities. The spectra are quite sim- 
ilar at all centralities: when fitting pT>5GeV/c with 
a power law function (oc p^), the exponents vary from 
n = -8.00 ± 0.12 in 0-5% to n = -8.06 ± 0.08 in the 80- 
92% (most peripheral) bin. Note that n = -8.22±0.09 in 
p+p collisions. The errors are combined statistical errors 
and systematic uncertainties. 

To quantify the comparison of spectra in heavy ion and 
p+p collisions, the nuclear modification factor (Raa) 



Raa = 



l/N cvt dN/dydp T 
(Tab) dapp/dydpT 



is used where a pp is the production cross section of the 
particle in p+p collisions, and (Tab) is the nuclear thick- 
ness function averaged over a range of impact parame- 
ters for the given centrality, calculated within a Glauber 
model 14]. Figure [2] shows Raa for ir° at different cen- 



tralities, the 0-5% bin is shown on Figure [3l The ref- 
erence p+p yield was obtained from the 2005 (Run-5) 
RHIC p+p measurement [18 1. 

Raa reaches ~0.2 in 0-10% centrality at pr>5GeV/c 
with very little (if any) j»t dependence. This trend is 
compatible with most current energy loss models but 
not with a semi-opaque medium assumption, where Raa 
would decrease with increasing j?t • While its magni- 
tude changes, the suppression pattern itself is remarkably 
similar at all centralities suggesting that the bulk Raa 
(integrated over the azimuthal angle) is sensitive only to 
the ./Vpart, but not to the specific geometry. Consequently, 
study of the px-integrated Raa vs. centrality is instruc- 
tive. 

Figure [4] shows the integrated nuclear modification fac- 
tor (pr>5GeV/c, and px>10GeV/c) for 7r°s as a func- 
tion of centrality, with the last two points indicating over- 
lapping 0-10% and 0-5% bins. In both cases the sup- 
pression increases monotonically with -/V par t without any 
sign of saturation, suggesting that larger colliding sys- 
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terns (such as U+U planned at RHIC) should exhibit 
even more suppression. 

The common power-law behavior (oc plj<) in p+p and 
Au+Au allows the suppression to be re-interpreted as 
a fractional energy loss S\ oss = 1 — Raa ^ wnere n 
is the power-law exponent, and we found that Si oas oc 
Fitting the integrated Raa with a func- 



tion i?AA = (1 - Sb^part) gives a = 0.58 ± 0.07 
for iV part >20 for p T >5 GeV/c, and a = 0.56 ± 0.10 
for p T >10GeV/c. The GLV @| and PQM [Hi mod- 
els predict that a w 2/3, which is consistent with the 
data. The fitted values of S are (8.3±3.3) x 10~ 3 and 
(9.2+4.9) xl0~ 3 for p T >5GeV/c and p T >10GeV/c, re- 
spectively. The fits are shown as dashed lines in Fig. |U 
Note that in this interpretation a constant S\ oss (indepen- 
dent of pt) implies that the energy loss increases withpx- 
We use the highest centrality (0-5%) -Raa data as 
shown on Fig. [3] to constrain the PQM model parame- 
ters. The procedure is described in detail in [171]. First 
we break up the errors of the measured points into Type 
A (pT-vmcorrelated, statistical © systematic, <jj), Type B 
(px-correlated, a bi , boxes on Fig. 2) and Type C (normal- 
ization, uniform fractional shift for all points, a c ). Then 
taking the theory curves calculated for different values 
of the input parameter p, one would normally perform a 
least-squares fit to the theory by finding the values of p, 
€b, e c that minimize: 



x 2 = 



E 



(Vi + tbVb, + tcViVc - Vi{p)) 2 



(1) 

where and e c are the fractions of the type B and C 
systematic uncertainties that all points are displaced to- 
gether. It is important to note that Eq. Q] follows the 
X 2 -distribution with n + 2 degrees of freedom when p, e b 
and e c are fixed, because it is the sum of n + 2 indepen- 
dent Gaussian distributed random variables. 

However, for the present data, the statistical and ran- 
dom systematic uncertainties are such that the shift in 
the measurement yt due to the correlated systematic un- 
certainties preserves the fractional type A uncertainty. 
Thus, we use a least squares fit of the quantity x 2 to 
estimate the best fit parameters, where \ 2 is Eq. [T]with 
o-j replaced by 5* = cr^j/j + e b a bi + e c yicr c )/y i: which is 
the uncertainty scaled by the multiplicative shift in yi 
such that the fractional uncertainty is unchanged un- 
der shifts. For any fixed values of e&, e c , x 2 still fol- 
lows the x 2 distribution with n + 2 degrees of freedom. 
The best fit, the minimum of x 2 , is found by standard 
methods (for example using a MINUIT type minimiza- 
tion algorithm) and the correlated uncertainties of the 
best fit parameters are estimated in the Gaussian ap- 
proximation by x 2 (eb,e c ,p) = Xmin + N 2 for N stan- 
dard deviation uncertainties. The right panel of Fig. 4 
shows the x 2 (et>, t c ,p) distribution for a wide range of 
values of the PQM model parameter (q) . Our data con- 



strain the PQM model transport coefficient (q) as 13.2 
^3 2 and ~t 5 ' 2 GeV 2 /fm at the one and two standard de- 
viation levels. These constraints include only the exper- 
imental uncertainties and do not account for the large 
model dependent differences in the quenching scenario 
and description of the medium. Extracting fundamental 
model-independent properties of the medium from the 
present data requires resolution of ambiguities and open 
questions in the models themselves which also will have 
to account simultaneously for the pt and centrality (av- 
erage pathlength) dependence. This work demonstrates 
the power of data for pion production in constraining the 
energy loss of partons. The data can be fitted with a con- 
stant in the entire pt> 5 GeV/c range as well: the slope 
of a simple linear fit is 0.0017 tomlt and ±0.0076 c/GeV 
at the one and two standard deviation levels. 

In summary, PHENIX has measured neutral pions in 
Au+Au collisions at ^/sjviv=200 GeV at mid rapidity in 
the transverse momentum range of l<px<20GeV/c, an- 
alyzing high statistics RHIC Run-4 run data. The shape 
of the spectra is similar for all centralities, as is the shape 
of .Raa(pt) at pT>5GeV/c. In central collisions the yield 
is suppressed by a factor of ~5 at 5GeV/c compared to 
the binary scaled p+p reference and the suppression pre- 
vails with little or no change up to 20GeV/c. Study- 
ing the integrated Raa vs. centrality we find that it 
does not saturate at this nuclear size; also the prediction 
Siogs oc TVp^jgl, [ll[ is consistent with our data and in 
this picture the energy loss increases with px- Using the 
0-5% (most central) Raa we find that the transport co- 
efficient (q) of the PQM model is constrained to 13.2 j^ 3 ' 2 
(-5 2) GeV 2 /fm at the one (two) a level. A simple linear 
fit with zero slope is also consistent with our data. 
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